> 

oo 



Domain wall of the totally asymmetric 
exclusion process without particle number 

conservation 



00 . 

^ ■ Yunxin Zhang*} $ 



Abstract 



O 

O 

t> . 

In this research, the totally asymmetric exclusion process without particle 
r S^ [ number conservation is discussed. Based on the mean field approximation and 

PL,- the Rankine-Hugoniot condition, the necessary and sufficient conditions of the 

■ existence of the domain wall have been obtained. Moreover, the properties of 

the domain wall, including the location and height, have been studied theoreti- 
cally. All the theoretical results are demonstrated by the computer simulations. 
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1 Introduction 

o 

J> . One-dimensional driven diffusion system is a very interesting research topic in 

^ recent years. They were shown to exhibit boundary induced phase transitions [T5] . 

and phase separation [9j [2TJ E] . In Ref . [17] , the effect of a single detachment site 
in the bulk of an asymmetric simple exclusion process (ASEP) was studied. In Refs 
[2T] [25] . the interplay of the simplest one- dimensional driven model, the totally asym- 
metric exclusion process (TASEP) with local absorption/desorption kinetics of single 
particles acting at all sites, termed "Langmuir kinetics" (LK) was considered. These 
models were inspired by the dynamics of motor proteins p] [2] p3] [27] [12] . which 
move along cytoskeletal filaments in a certain preferred direction while detachment 
and attachment can also occur between the cytoplasm and the filament, and, in a 
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very different setting, by dynamics of limit orders in a stock exchange market. Being 
an equilibrium process, LK is well understood, while the combined process of TASEP 
and LK showed the new feature of a localized domain wall in the density profile of 
the stationary state [21]. 

The TASEP is defined on a one-dimensional lattice of size N. Each site can either 
be empty or occupied by one particle. In the bulk, particles can hop from site % to site 
2 + 1 with unit rate, provided the target site is empty. At site 1, particles can enter the 
lattice from a reservoir with density a , provided the site is empty. They can leave the 
system from site N into a reservoir of density (3 with rate 1 — (3 . Thus in the interior 
of the lattice, the particle number is a conserved quantity. The phase diagram and 
steady states of the TASEP as a function of the boundary rates are known exactly 
[T6] [21] [6]. Furthermore, a theory of boundary induced phase transitions exists, 
which explains the phase diagram quantitatively in terms of the dynamics of domain 
wall [14J. In the field of TASEP, Joel L. Lebowitz et al have done many excellent 
research and obtained much useful theoretical results [El E] • 

Similar as in [25] , in this research we equip the system with the additional feature 
of local particle creation at empty sites with rate uj a and annihilation with rate Ud- 
In the thermodynamic limit iV — ► oo, the case of the local rates being of the order 
of 1/N is the most interesting one |22j. It turns out that the presence of the kinetic 
rates significantly change the picture of TASEP, producing a completely reorganized 
phase diagram. In [2T], the authors showed by computer simulations and mean-field 
arguments that, in this nonconserved dynamics, one can have phase coexistence where 
low and high density phases are separated by stable discontinuities (domain wall) in 
the density profile. Recently, this dynamics was also studied theoretically in [18] [19]. 

Up to now, the properties of the totally asymmetric exclusion process without 
particle number conservation have not been well studied theoretically, though many 
of which have been found in the computer simulations. It is no doubt that some 
particular properties of this process are very difficult to be found only by computer 
simulations, and the right theoretical analysis should be given to the full understand- 
ing of which. In this research, the necessary and sufficient conditions of the existence 
of domain wall and the properties of the domain wall, for this totally asymmetric ex- 
clusion process without particle number conservation, will be theoretically discussed. 
Two basic questions will be answered theoretically: (1) when and where does the 
domain wall exist? (2) How do the location and height of the domain wall change as 
the parameters change? 

This paper is organized as following. In the next section, the mathematical model 
and some basic results of this process are introduced. The necessary and sufficient 
conditions of the existence of the domain wall will be given in section 3, then the 



2 



properties of the domain wall, including the location and the height, will be discussed 
theoretically in section 4. Finally, some concluding remarks are given in the last 
section. 



2 Mean field approximation 

The equations of the bulk dynamics of the totally asymmetric exclusive process 
with particle creation and annihilation are the following P, [21] 



din ' 

= 7li_i(l - 7li) - rii(l - n i+ x) + u a (l - Hi) - u d rii 



2 < i < N 



while at the boundaries 



~dt ' 

driN 

~dT 



a(l - m) - m(l - n 2 ) 
-- n N -i(l - n N ) - (3n N 



(2) 



where the occupation numbers = 1 for a site occupied by a particle and = 
for an empty site. For fixed total length L = 1 and N — ► oo one gets the differential 
equation for the average profile in the stationary state 



(3) 



(2 P (x) - i)p'(x) = (n a + n d ) P (x) -n a o<x<i 



where the reduced rates Q a = Nu a , Q d = Nu d (because of the particle- hole symmetry, 
we restrict the discussion to the case fl a > fl d [21]). In the following, let pi a and p T p 
(or pi and p r for simplicity) be the solutions of Eq. (J3j) with boundary conditions 
Pza(O) = a and p r p{l) = 1-/3 respectively. It is to say, p\ a and p r p satisfy the 
following two equations respectively (see [TP] ) 



2{p la - a) K-l 



In 



K-(K + l)p 



la 



K - (K + l)a 



x 



+ 



K 



In 



K-{K + !){!- (3) 



K-(K + l)p 



r/3 



1 — X 



(4) 
(5) 



where K :— > 1. 



3 The existence of the domain wall 

Due to the Rankine-Hugoniot condition |4], at the location x$ of the domain wall 
ofEqn. (ED, 

(^-p + ) = (p 2 (6) 
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Figure 1: For < a < 0.5 and 0.5 < 1 - (3 < if p£{fl) < l,p~p(l - a) < 0, 

there exists domain wall in (0,1) (left). For < a < 0.5, < 1 — /3 < 1: if 
Pia(P) < 1 an d ~ a ) — 0; there exists domain wall in (0, 1) (right). 

should be satisfied, where p_ = lim p(x), p + = lim p(x). It can be easily found 
that the condition can be simplified as 

P- + P+ = l (7) 

Therefore, to the Eqn. ([3]) with boundary conditions p(0) = a and p(l) = 1-/3, 
there exists domain wall in the interval (0, 1) if and only if there exists a location 
< Xs < 1 at which p_ + p + = 1, moreover, 



p(x) 



pi(x) < x < xs 

p r {x) Xs < X < 1 



From the conditions (jSJ) and the Eqns. dl]) ©, it is not difficult to obtain the 
following theoretical results: 

(I) For < a < 0.5, < 1 - p < the necessary and sufficient conditions of the 
existence of domain wall are (see Fig. [1] (left)) 

Pte(7)<l and p^ x (l-a)<0 with 7 = min(0.5, p) (9) 



(II) For < a < 0.5, < 1 - p < 1, the necessary and sufficient conditions of the 
existence of domain wall are (see Fig. [TJ (right)) 

PTM<1 and p^(l-«)>0 (10) 

(III) For 0.5 < a < 1, there no domain wall exists in (0, 1) (for details, see [26J). 
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Figure 2: The relationship between the location of the domain wall and the parameter 
Qd, for < a < 0.5, -^-j- < (3 < 0.5 (left). xs(fld) is monotonously decreased as a 
function of Q d for < a < 0.5, ^ < (3 < 0.5 (right). 

4 The properties of the domain wall 

In this section, we will discuss the properties of the location x$ and the height 2e 
of the domain wall, which can be regarded as functions of the parameters K, fid, a, (3. 
Where 2e = \p + — p_| = \2p + — 1|. In the following, we assume < a < 0.5 (which is 
the necessary condition of the existence of domain wall), < (3 < 0.5 (if /3 > 0.5 and 
the domain wall exists, it is equivalent to the case in which (3 = 0.5, see [26J), Qd > 
and K > 3. For the sake of the convenience, we define the following functions 

A = A(e, K) := 2e + 2Ke — K + 1 B = B(e, K) := 2e + 2Ke + K - 1 

C = C(K 1 (3) := K(3 + (3 -1 D = D(K, a) := Ka + a — K (11) 

E = E(e, a) := 1 - 2a - 2e F — F(e, (3) := 1 - 2(3 - 2e 

4.1 The Properties of the Location of the domain wall 

From Eqns. (jlj) and ([5]), we can get the following theoretical results: 

(a) 

dxs_ = A(e,K)-Aex s (K + l) = 2e(K +!)(! — 2x s ) — (K — l) 

dn d 4ettd(K + l) Aen d (K + l) [ ' 

For ^ < (3 < 0.5, i.e. 0.5 < 1 - (3 < it can be easily proved that 2e < f^, 
which implies A(e, K) < 0, so |^ < 0. It is to say that the location x$ of the domain 
wall is monotonously decreased as a function of the parameter fid (the corresponding 
computer simulations are plotted in Figure[2]). For < (3 < -j^tt, i.e. jM-r < l—(3 < 1, 



K+1 , l-^- K+1 

dx s = 2e(K+l)(l-2x s )-(K-l) 
dfl d AeVl d (K+l) 



(13) 
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Figure 3: The relationship between the location of the domain wall and the parameter 
Q d : < a < 0.5, < (3 < (left). x s (Q d ) is monotonously decreased as a function 
of fid for < a < 0.5,0 < {3 < -^-j (right). In all the computer simulations, 
2e(K + 1)(1 - 20:5) < (K - 1) is satisfied. 

if 2e(K + 1)(1 — 2xs) < {K — 1), the location xs of the domain wall is also decreased 
as the parameter fl d increases (Figure [3]). Otherwise, the location xs of the domain 
wall is increased as the parameter Q d increases (Figure HI), 
(b) 



dxs 1 



dK A{K+l)e 



h ~ 3 {(A + B)x s -A} 



K 2 -V 



\{K 2 - l)(2a - 1) + Ae(K + 1)D]CE - [Ae(K + l)C + 2(K + 1)(1 - 2(3)]DF~ 



(14) 



tt d (K+ l) 2 (K- 1)CD 

For < a < 0.5, < (3 < 0.5, it can be verified that 

A<0 A + B>0 C>0 D<0 E>0 F<0 (15) 

which imply 

^ < for K > 3 (16) 

i.e. the location xs of the domain wall is monotonously decreased as a function of 
the parameter K for K > 3, < a < 0.5, < [3 < 0.5 (Figure [5]). At the same 
time, Eq. ([TJJ can be reformulated as 

W = - 4(K\r )f n d (K + mK-i)CD «* - ^<- K + " J > + & - « 

+ 8e(/3 - a) + (K 2 - l)(2a - 1)CE + (if - l)(2/3 - 

(17) 

For < (3 < jA^, it is easy to verify 

C < D < £ > F > (18) 
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Figure 4: The relationship between the location of the domain wall and the parameter 
&d for < a < 0.5,0 < f3 < jA^' monotonously increased as a function of the 
parameter when Qd is small, then monotonously decreased when is large enough 
(left). The corresponding figure of the function xs(^d)- (right). At the critical 
point, 2e(K + 1)(1 - 2x s ) = (K - 1). 




Figure 5: For < a < 0.5, -^rj < (3 < 0.5 the location xs of the domain wall is 
decreased as the parameter K increases (left). The figure of the function xs{K) for 
< a < 0.5, < 13 < 0.5 (right). 
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Figure 6: For < a < 0.5, < (3 < 2e(K + 1)(1 - 2x s ) < (K - 1) the location 
of the domain wall is decreased as the parameter K increase (left). The figure 
of the function x s {K) for < a < 0.5, < (3 < 2e(K + 1)(1 - 2a; 5 ) < (K - 1) 
(right). 

if 

2 £( * +1 )(i-2. S ) < ( k-i )+ 8 ^ - °> + ^ - ^ 2a - i£g + g - agg - 1 > gF 

i\ — 3 

(19) 

then |^ < 0, the location xg of the domain wall is decreased as the parameter 
K(> 3) increases (Figure ED- Since 

(K 2 -l)(2a-l)CE + (K -1)(2(3-1)DF 1 

* ^ — — — ^° °^'^Ffi (20) 

the breakdown of the inequality (TT9~j) is difficult to be found in the computer simula- 
tions, 
(c) 

^ ^i'- 2 °!n ^ *» V0 < a, 13 < 0.5 (21) 
da AeVt d (K + 1)D v ' 

so the location 25 of the domain wall is monotonously decreased as a function of the 

parameter a for V K > 1,0 < a < 0.5,0 < (3 < 0.5 (the results of the computer 

simulations are plotted in Figure [3, Figure [8]) 

(d) 

in fact, A < Q,C > ii < {3 < 0.5; and A > 0, C < if ^ < /3 < 0.5. So 
the location of the domain wall is monotonously increased as a function of the 
parameter f3 for V K > 1, < a < 0.5, < (3 < 0.5 ( Figure El Figure M). 



8 





Figure 8: The location xs of the domain wall is decreased as the increase of the 
parameter a for < (3 < (left). The figure of the function xs(cx) for < (3 < -^-j- 
( right). 




Figure 9: The location x§ of the domain wall is increased as the increase of the 
parameter (3 for -^-j < (3 < 0.5 (left). The figure of the function xg{j3) for -^-^ < 
P < 0.5 (right). 
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Figure 10: The location xs of the domain wall is increased as the increase of the 
parameter (3 for < (3 < -^q-j- (left). The figure of the function xs{(3) for < (3 < -^-j- 
(right). 

4.2 The Properties of the Height of the domain wall 

In view of the Eqns. (JIJ) ([51) and the definition of the height of the domain wall 

2e, using the chain rule of the derivative, we can get the following theoretical results: 

(a) i 

., IP | >0 for — <(3<0.5 

J - K + l- P - (23) 



dn d 16{K + l)e> <0 for 0<(3< 



K+1 

It can be verified that A < 0, B > for < (3 < 0.5; and A > 0, B > for 
< f3 < ^Vj- So the height 2e of the domain wall is monotonously increased as a 
function of the parameter Q d for -^q-j- < /3 < 0.5, and monotonously decreased as 
a function of the parameter Q d for < /3 < (Figure [2] (left), Figured (left), 
Figure H (left) , Figure [TTJ . 
(b) 

9e f [(/<" - l) 2 + 2BD)ACE + [(^ - l) 2 + 2AC]BDF 

dK ~~ { (K+1) 2 (K -1)CD 

(K-3)ABtt d \ 1 



(24) 



+ K 2 - 1 J 16(if + l)e 2 
where 

(K - l) 2 + 2BD = (K 2 - l)(2a - 1) + 4e(ir + 1)D 
(K - l) 2 + 2AC = Ae{K + l)C + (K - 1)(1 - 2/3) 
From (15) and (24), one can know 



(25) 



<9e 1 

> for K > 3, — — - < (3 < 0.5, < a < 0.5 (26) 



_ — K + 1 
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Figure 11: The figure of the function 2e(f2d) for -^q-j- < (3 < 0.5 (left). The figure of 
the function 2e{Q d ) for < (3 < (right). 




Figure 12: The figure of the function 2e(K) for ^ < (3 < 0.5 (left). The height 
2e of the domain wall is not a monotone function of the parameter K for 1 < K < 
3, < (3 < j^, 0<a< 0.5(right). 



so the height 2e of the domain wall is monotonously increased as a function of the 
parameter K for K > 3, j^— < (3 < 0.5, < a < 0.5 (Figure [5] (left), Figure 
T2l (left)). However, for < /? < ^Vj, the height 2e of the domain wall is not a 
monotone function of parameter K (Figure [T2l (right) ) . 
(c) 

de (1 - 2a)AB [ ^ ° F+~T - ^ - °' 5 

da 16(Jjf + l)VD] <Q Q< 5< 1 

I ~ ~ M ~ K + l 

so the height 2e of the domain wall is monotonously increased as a function of the 
parameter a for ^ < /3 < 0.5, and monotonously decreased as a function of the 
parameter a for < (3 < -^q-j- (Figure [7|( left), Figure [8] (left), Figure [TBI. 



(27) 




0.16 0.18 0.2 0.22 0.24 0.26 0.28 0.3 0.32 



3 

I 0.75 



K=3, a d =0.08, 


=0 




12^=11. Uo. p=U. I 


K=3, £l d =C 


.08, p=0.1 








K=2, £1 =0.15, p=0- 


k=3, n d = 


0.15, p=0.1 —~ 



Figure 13: The figure of the function 2e(a) for -^-j- < (3 < 0.5 (left). The figure of 
the function 2e(a) for < (3 < -^-j- (right). 





K-2, £1 =0-5, a=0.05 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 



Figure 14: The figure of the function 2e(/3) for -^-j- < (3 < 0.5 (left). The figure of 
the function 2e(J3) for < (3 < ^_ (right). 



(d) 



de (1 - 2(3)AB 



< for < a, /? < 0.5 



(2f 



9/3 16(K+l)e 2 C 

since ^ < and i? > for V < a, /3 < 0.5. Therefore, the height 2e of the domain 
wall is monotonously decreased as a function of the f3 for V0 < (3 < 0.5 (Figure [9] 
(left) , Figure M (left) , Figure ED . 



5 Conclusions 

In this research, the totally asymmetric exclusion process without particle num- 
ber conservation in large particle number limit have been studied theoretically. Two 
questions are answered completely: (1) when and where does the domain wall exist? 
(2) How do the location and height of the domain wall change as the parameters 
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a, (3, K, fid change? (see Table 1, where j (j) means the function is a monotonously 
increased (decreased) one, || (if) means the function has an unique maximum (min- 
imum) point, |? means that the function might be monotonously decreased, but the 
proof is not completed in this research) 

In summery, we have found: (1) For < a < 0.5, < 1 — (3 < j^tj, there exists 
domain wall if and only if p^(l) < 1 and Pr^O- ~ a ) < 0, where 7 = min(0.5, (3). 
(2) For < a < 0.5, < 1 — (3 < 1, there exists domain wall if and only if 
Pi 1 {ft) < 1 an d P^ 1 (l — «) > 0. (3) For 0.5 < a < 1, the domain wall doesn't ex- 
ist. (4) The location xg of the domain wall is monotonously increased (decreased) as 
a function of the parameter f3 (a). The height 2e of the domain wall is monotonously 
decreased as a function of the parameter (3. (5) For < (3 < 0.5, K > 3, The 
location x$ of the domain wall is monotonously decreased as a function of the pa- 
rameters fid and K; the height 2e of the domain wall is monotonously increased as a 
function of the parameters fid and K. (6) For < (3 < ^tj, , the height 2e of the 
domain wall is monotonously decreased as a function of the parameters fid and a. 



< a < 0.5 




if (> 3) t 


a | 


/ST 


7^T < P < 0-5 


^5 1, 2e T 


x S i, 2e T 


xs i, 2e T 


^5 T, 2e I 


0</?< xTT 


x 5 TI, 2e | 


x 5 I?, 2e IT 


x 5 I, 2e I 


^5 T, 2e | 



Table 1: The properties of the domain wall. 



Recently, totally asymmetric exclusion processes with internal states and particle 
detachment and attachment, which serve as generic transport models in various con- 
text, have been introduced and extensively studied ([23] [11] [20J). Using the similar 
methods as in this research, these generalized models also can be studied theoretically, 
and the corresponding results will be given in the future. 
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